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We present a lattice model to study the equilibrium phase diagram of ordered alloys with one mag- 
netic component that exhibits a low temperature phase separation between paramagnetic and fer- 
romagnetic phases. The model is constructed from the experimental facts observed in Cua-^AlMna, 
and it includes coupling between configurational and magnetic degrees of freedom which are ap- 
propriated for reproducing the low temperature miscibility gap. The essential ingredient for the 
occurrence of such a coexistence region is the development of ferromagnetic order induced by the 
long-range atomic order of the magnetic component. A comparative study of both mean-field and 
Monte Carlo solutions is presented. Moreover, the model may enable the study of the structure of 
the ferromagnetic domains embedded in the non-magnetic matrix. This is relevant in relation to 
phenomena such as magnetoresistance and paramagnetism. 



I. INTRODUCTION 

Recently, renewed interest has been addressed to fer- 
romagnetic ordered alloys. This is because of the 
unique properties arising from the interplay between 
elasticity, magnetism and (configurational) atomic or- 
der. From the point of view of applications, the de- 
velopment of new actuator materials having very large 
magnetostrainsa is particularly interesting . Also-relevant 
is the possibility of haying superparamagnetisrrjj and gi- 
ant magnetoresistanceB both associated with coexistence 
of magnetic domains (large magnetic particles) embed- 
ded in a non-magnetic matrix. This mixed phase has 
been observed, for instance, in the Cus-^AlMn^ Heusler 
alloy. 

The Heusler alloys are ternary intermetallic com- 
pounds with the composition X2YZ and a low temper- 
ature L2\ structure. At high temperatures the sta- 
ble phase corresponds to a disordered bcc lattice, also 
called A 2 phase, and undergoes a two-stage disorder- 
order transition A2 — > B2 — ► L2%, as the temperature 
is decreased. -Especially interesting are the Mn-based 
Heusler alloysn 0, which exhibit a magnetic moment ap- 
proximately located on the Mn atomsH. Amonffpthem, 
the most extensively studied are the Ni 2 GaMnBlia and 
the Cu2AlMnErtZl alloys. In both cases, the L2\ phase 
is ferromagnetic but the B2 phase is paramagnetic. This 
close relation between atomic order and magnetic-proper- 
ties has been known to scientists for many yearala. Addi- 
tionally, these alloys exhibit shape-memory effects, inti- 
mately related to the structural transition, of the marten- 
sitic typeO, undergone at low temperatures. It has been 
suggested that the control of shape-memory properties 
by application of an external magnetic field is a principle 
for operation of the new class of actuator materials^'ElL 

In Cu-Al-Mn the martensitic transition only existstil 
for compositions which are very far from the stoi- 
chiometry (Cu^AJMn) where the L2i ordered phase 
is paramagneticcd. Nevertheless, the influence of 
magnetism coming from Mn is revealed in several 



experimentsE3. As well as the phase transitions men- 
tioned above, the system exhibits, at low tempera- 
tures, a spinndal decomposition along the line Cu 3 Al - 
Cu2AlMnOli3cil. We will centre our attention on this 
two-phase region and denote the Cu-rich portion of the 
phase diagram of interest in this paper by Cus^xAlMn^, 
with < x < 1, . In Figure |l| we show schematically the 
corresponding phase diagram as it is obtained from ex- 
periment. The continuous lines are drawn from the data 
in Ref. 24, whereas the points at x — and x = 0.28 are 
from Refs. |2fj and |26| respectively. 

The low temperature ordered structures for the lim- 
iting values of x are different- The CU3AI binary alloy 
is DO3 at low temperaturesj23 whereas the Cu^AlMn is 
L2i and ferromagnetic, with a relatively high Curie tem- 
perature (~ 630 K)ES The ferromagnetism of the L2i 
phase appears as a consequence of the atomic ordering 
of the Mn atoms. In this sense it is known that prop- 
erties such as the saturation magnetic moment depends 
on the degree of order of the Mn atomsca It then natu- 
rally follows that the absence of magnetism (long-range 
magnetic order) either in the high temperature i?2-phase 
or in the low temperature phase (DO3 or L2\), for small 
values of x, might well be related to the tendency for 
the Mn atoms to distribute themselves randomly at the 
different lattice sites. On the other hand, by increasing 
the amount of Mn, for instance in the Cu3AlMn2 allpM, 
the resulting magnetic interaction is antiferromagneticEa 
Such different magnetic behaviour may be understood in 
terms of the oscillatory RKKY intepictipH between the 
magnetic moments of the Mn atomsc2lllIle2l. 

The phase separation or miscibility gap in Cu2AlMn 
occurs at temperatures below ~ 600KE3 (see Figure |l|) 
and gives rise to a coexistence region between a non- 
magnetic phase and a ferromagnetic phase, with low 
and high Mn conteBt respectively. The occurrence of 
superparamagnetismo or magnetoresistancea are both di- 
rectly related to the existence of magnetic clusters ( sta- 
ble L2i domains) immersed in the non-magnetic (DO3) 
matrix. Some aspects of this phase diagram are not to- 



1 



tally clear. Firstly, the persistence of a stable DO3 phase 
for small values of x. Kainuma et ai.O, by using X- 
ray diffraction measurements, have detected an abrupt 
change in the intensity of the superstructure peaks at 
x ~ 0.32. It should be mentioned that this effect was not 
found in other earlier studiesliil. Other important infor- 
mation, not yet available, refers to the different atomic 
distributions for the non-stoichiometric L2\ structure. 
Some assumptions on this matter will be required in or- 
der to perform a theoretical study. Other aspects that 
need to be discussed refer to the characteristics of the co- 
existing phases. They will depend on the location of the 
DO3-L21 (according to the results in Ref. ^) and mag- 
netic transition lines with respect to the coexistence line. 
More precisely, depending on the temperatures at which 
such interphases end on the coexistence line, the phases 
may be different in atomic order (D03,L2\) or/and mag- 
netic order (ferromagnetic, paramagnetic). In this sense, 
even the coexistence of two different paramagnetic L2i 
and L2[ phases (upper part of the miscibility gap in Fig.- 
ure 1), with a very similar content of Mn, are suggesteco. 

In this paper we present a lattice model able to repro- 
duce the main features of the equilibrium phase diagram 
in this two-phase region. The details of the model will 
be derived from a microscopic description of the atomic 
and magnetic properties of Cu3_ a AlMn x alloys. Nev- 
ertheless, it can be applied to other systems. Practi- 
cal reasons will require several hypotheses which in some 
cases are not totally justified a priori but only later from 
agreement of obtained results with experimental data. 
This agreement is indicative that the model captures 
the essential physics and provides a starting point for 
future more exhaustive studies. The model is a projec- 
tion of the ternary alloy onto a binary system, when one 
of the species is magnetic. It is constructed on the ba- 
sis that the main physics of the phenomena lies on the 
atomic ordering of the magnetic component which more- 
over is taken to be always the less abundant. The ef- 
fective Hamiltonian accounts for a purely configurational 
ordering energy between first neighbouring pairs so that 
at low temperatures the magnetic atoms tend to be sec- 
ond neighbours. Then a simple ferromagnetic pair in- 
teraction between next-nearest neighbours is enough to 
give rise to a low temperature phase separation between 
a non-magnetic phase and a ferromagnetic phase that, 
moreover, may have different ordered structures. 

It has been suggestcdcil that the occurrence of the two- 
phase region in Cua-^AlMnj, cannot be attributed to ei- 
ther chemical(configurational) or magnetic ordering. In 
this work, we use a very simple microscopic model to 
show that the coupling between both atomic (configura- 
tional) and magnetic orderings is sufficient to give rise 
to a decomposition between two phases at low temper- 
atures. This coupling operates in such a way that as 
the atomic ordering develops the (indirect) exchange in- 
teractions between the atomic moments of the magnetic 
particles produce a long-range ferromagnetic order. 

The remainder of the paper is organized as follows. In 



section II we introduce the model. Section III is devoted 
to its mean-field solution. In order to better understand 
the nature of the different phases and the behaviour of 
several measurable quantities we also solve the model by 
using Monte Carlo numerical simulations. This is pre- 
sented in section IV. Finally in section V we summarize 
and conclude. 



II. MODEL 

In the present study, our main goal is to understand the 
formation of the miscibility gap in Cu3_ x AlMn x along 
the line < x < 1. The complexity inherent to the de- 
scription of a magnetic ternary alloy has led us to make 
simplifications that we shall discuss in this section. In- 
deed, the quest for reasonable simplifications becomes 
compulsory in order to perform the Monte Carlo numer- 
ical simulations. Although the inclusion of too many in- 
gredients (and thus free parameters) in the model may 
lead to a better fit of the available data (in our case 
scarce), it may hide the understanding of the relevant 
physical mechanism underlying the phase diagram prop- 
erties, which we believe is the coupling between the long- 
range configurational (chemical) ordering and the mag- 
netism of the Mn atoms. 

The equilibrium structure of Cua-^AlMnj, can be de- 
scribed as an underlying bee structure formed by the 
superposition of four interpenetrated fee sublattices, 
named a, (3, 7 and 5 (see Figure ||a). In order to de- 
scribe the different phases of the system it is convenient 
to specify the occupation probabilities p x of the different 
species X (=Cu,Al,Mn) in the four different sublattices 
£ (=a, (3, 7, 8). Tables | and |l| summarize the occupa- 
tion probabilities for the limiting DO3 (x — 0) and L2i 
(x = I) stoichiometric phases. For intermediate values of 
x a more elaborate discussion is required. 

We start with the region corresponding to small values 

of x, x > 0. Recentlyu, X-ray diffraction experiments 
shown that the DO3 structure persists above the coexis- 
tence region for values of x up to 0.32. In other words, 
the addition of a small amount of Mn does not break the 
symmetry p x =p x = Px ■ This is an (a priori) unexpected 
result, given the different atomic environments of these 
sublattices in the DO3 phase. In any case, it seems clear 
that entropy plays a very important role in the stability 
of this homogeneous DO3 phase. From Figure [l] it fol- 
lows that for x > the stability is extended to higher 
temperatures as the value of x increases. A natural hy- 
pothesis is, therefore, to assume that (for low values of 
x) the Mn atoms behave as impurities that are randomly 
distributed on the four different sublattices. The cor- 
responding occupation probabilities p x are indicated in 
Table |n| 

In the x < 1 region the stable phase is of the L2i 
type. There are several atomic distributions that are 
compatible with the symmetry p x —p x ^ p x p x . Ta- 
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blc IV displays the occupation probabilities in the most 



straightforward case for which the Mn concentrates in a 
unique sublattice. Alternatively, in a more general way, 
one might write the occupation probabilities (Table |V|) 
in terms of a free parameter A (0 < A < 1). Notice 
that these atomic distributions account for a continuous 
change from the DO3 phase (A = 0) to the L2\ phase 
(A > 0). 

The next step is to introduce the two major simplifi- 
cations of the model: 

1. The structures described in tables I to V have the 
existence of two n.n. sublattices (7 and 5) in com- 
mon which contain most of the Cu atoms and have 
identical occupation probabilities: p x = p x , \/X. 
Experimentally, this symmetry with respect to the 
7 and S sublattices seems to be satisfied for any 
concentration and temperature ranges. From now 
on we forget about them and concentrate on the 
atomic distribution behaviour on the other two re- 
maining sublattices, motivated by the feeling that 
the breaking down of the a-(3 symmetry is crucial 
in the ordering of the Mn atoms at low tempera- 
tures. This, of course, will restrict the validity of 
our study to temperatures below the B2-DO3 tran- 
sition, which is precisely the region of interest here. 
Therefore, the model will be defined on a simple 
cubic lattice divided into two sublattices, a and /3, 
as illustrated in Figure §b. 

2. Continuing with our assumption that the main 
physics lies on the atomic ordering of the Mn 
atoms, we shall proceed further by distinguishing 
between magnetic and non-magnetic atoms only. In 
our binary alloy model, Ai_ c B c , the non-magnetic 
species A stands either for Cu or Al, whereas the 
magnetic species B stands for Mn and the composi- 
tion is restricted to c < 0.50. The behaviour of the 
B atoms on sublattices a and (3 can be regarded as 
a simple order-disorder transition. For small val- 
ues of c both sublattices are equally populated by 
B-atoms (behaving as impurities) while for larger 
values of c, B atoms occupy preferably one of the 
two sublattices. Moreover, this behaviour depends 
on temperature. As regards the configurational or- 
dering, the model gives rise to two phases only: dis- 
ordered and ordered, corresponding to low (DO3) 
and high (L2i) content of the magnetic species re- 
spectively. Keeping this correspondence in mind, 
in what follows we shall use the simplified notation 
D (disordered) and O (ordered). 

We notice that the quantitative study of properties 
such as the magnetization, susceptibility or other prop- 
erties related to magnetism (magnetoresistance, etc..) is 
not our goal here. Rather, we shall focus on how the de- 
velopment of long-range ferromagnetic order (resulting 
from the interplay with the atomic order) determines the 
phase diagram at low temperatures as a function of the 



content of the magnetic species c. In what the model 
description concerns, this can be achieved by considering 
localized Ising-like spin variables Sj = ±1 associated with 
each B atom. 

We start from the following pair-interaction effective 
Hamiltonian: 



n = n c + n m = 

w 

= Y.[ N AAtAA + N k AB e k AB + X< i:ll < 



(1) 



fc=l 



fc=i 

where TL C and TL m are the configurational and magnetic 
energy contributions respectively. The summation is 
performed over the different fc-nearest-neighbour shells 
(up to k = w), N xx is the number of fc-th nearest- 
neighbour X-X pairs and e xx are their corresponding 
pair-interaction energy. Note that the magnetic contri- 
bution only involves B atoms. We have indicated by B + 
and B~ the two possible magnetic states. To preserve 
the symmetry under exchange of the + and — magnetic 
states we take e+ + = e k __ Vfc. 

Following standard procedures, we write Hamiltonian 
(|l|) in terms of Ising-like variables defined at each lattice 
site. Let us index the sites of the cubic lattice by i — 
1, • • • , N (N = L x L x L). At each lattice site i we 
define the following two coupled two-state variables Oi 
and Si . The variable Cj = +1,-1 represents the non- 
magnetic and magnetic species (A and B) respectively, 
then, provided <Ji = —1, we define Si = +1,-1 describing 
the two possible magnetic states of each B atom. 

Considering interactions up to next-nearest neighbours 
(w = 2), the configurational energy term in equation (Q) 
can be written, neglecting constant terms, as: 

nn nnn N 

H c = J{ E + J 2 C £ aw +E C Y, °i > (2) 

ij ij i=l 

where the first two summations are extended to near- 
est neighbour (n.n.) and next-nearest neighbours (n.n.n.) 
respectively, and the Hamiltonian parameters are: 
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(3) 
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where z\ = 6 and z 2 = 12 are the number of n.n. and 
n.n.n. of each lattice site respectively. In the Canonical 
ensemble, the last term in equation (0) is just a simple 
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energy shift which depends on the alloy concentration. 
As regards the magnetic energy term it can be rewritten 
as: 
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(6) 
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(7) 



(8) 



We notice that the latter two terms in equation (||) do 
not depend on the magnetic variables {Si}. Expanding 
the different contributions in (^|) and ignoring constant 
terms, the Hamiltonian becomes: 



H = ( Jf + K?) + (J2 + K?) °W 



(9) 



n.n.n. 



The superscripts in the model parameters denote its con- 
figurational ( c ) or magnetic ( m ) origin, whereas the sub- 
scripts mean first- (i) or second- (2) neighbour interac- 
tions. In order to reduce the number of free model pa- 
rameters we set J™ = 0. Indeed, the n.n magnetic in- 
teraction between B — B pairs is not essential for our 
present purposes since we restrict ourselves to the case 
in which the ferromagnetism develops in the configura- 
tionally ordered phase. Furthermore, by using reduced 
energy units H* = H/(Jf + K™), we get the following 
minimal model Hamiltonian: 



(10) 



K* = - 



Jf + K™ ' 



(11) 



which measures the ordering energy between second- 
neighbour pairs either A — A, B — B and A — B, in- 
dependently of the magnetic state of atom B, and 



T* 
Jo 



Jf + K{ n 



(12) 



which accounts for the ferromagnetic interaction between 
second-neighbour B — B pairs. 



III. MEAN FIELD SOLUTION 

This section is devoted to the solution of the model 
introduced previously for the Ai_ c B c binary alloy by us- 
ing standard mean-field techniques based on the Bragg- 
Williams approximation. We denote the occupation 
numbers for each component (X = A, B + , B~) in each 
sublattice (E = a, (3) by and consider the following 
order parameters: 
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(13) 



(14) 



(15) 



where c (0 < c < 0.5) is the molar fraction of the mag- 
netic species, 77 (0 < 77 < 2c) is the atomic order param- 
eter and m (0 < m < 2c) measures the magnetization 
of the system. Using standard procedures, in the Grand 
Canonical formulation, we obtain the following expres- 
sion for the internal energy: 



E = NJ Zl 



\{l-2cfll-K* Z -l]~\ll + K* Z -l] 
2 z\ 2 z\ 



(16) 



/i*(l-2c)- jf— m 2 

Zl 



where J = Jf + K™ and fi* is the chemical potential 
difference between the two species. The corresponding 
entropy is given by: 
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Expressions ( |l6| ) and (17) produce the following free en- 
ergy: 



where the parameters are: 
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+ (c + | + 2m) ln(c + | + 2m) 
+ (c + | - 2m) ln(c + - - 2m) 
+ 2(c-|)ln(c-|)-4cln2 

with T* = 2^. and J > 0. The free energy in @, in 
the absence of magnetism, reduces to the standard case of 
order-disorder, but one of the species is twice degenerate. 
When magnetism is taken into account, model (|l^) ex- 
hibits two phase transitions respectively associated with 
the order parameters 7/ and m. We denote the respective 
transition temperatures by T* and T^. Since we are in- 
terested in the case of T* n < T* the model parameters 
must be taken so that (1 + K *) > J| > 0. 

The equilibrium temperature dependence of the order 
parameters was obtained from direct minimization of the 
function ( |l8| ) . In figure |^ we show the fi* — T* section of 
the phase diagram for different values of JJ = 0.60, 0.75 
and K* = 0.0, 0.20 and 0.40. Three different phases may 
appear. The Disordered-Paramagnetic (DP) phase with 
7] = and m = 0, the Ordered-Paramagnetic (OP) phase 
with 77^0 and m = and the Ordered-Ferromagnetic 
(OF) phase with 77 ^ and ra/0. Both parameters Jf 
and K* have the effect of increasing the stability of the 
ordered (OP and OF) phases. Continuous lines stand 
for second-order phase transitions, whereas the dashed 
ones stand for discontinuous phase transitions. The in- 
tersection between the three interphases corresponds to 
a bicritical point in cases (a) and (b) , whereas for (c) and 
(d) it corresponds to a triple point. The DP-OF transi- 
tion is always first order, whereas the other two OP-OF 
and DP-OP may be second or first order. When the tran- 
sition is first order, a phase separation shows up in the 
c — T* section. This is illustrated in Figure [| for cases 
b) and (c) corresponding to the previous picture (Figure 



In both cases of Figure |] a phase separation be- 
tween a non-magnetic (paramagnetic) and a ferromag- 
netic phases exists. At low temperatures the coexist- 
ing phases (DP+OF) are also different in their atomic 
ordered structure, whereas at moderate temperatures 
(OP+OF) both exhibit the same atomic structure. Be- 
sides, for case (b) (a larger value of K*) a phase sepa- 
ration (DP+OP) between two non-magnetic phases ap- 
pears. In this case, there exists a line of triple points 
(horizontal dot-dashed line). In Figure ^ we show the 



corresponding temperature behaviour of the order pa- 
rameters 77 and to for different values of the composition 
c. This information is obtained from the calculations 
presented in Figure || taking into account the fact that in 
the phase separation region the system is heterogeneous 
and that at constant concentration both the character- 
istics and the amount of the coexisting phases change 
with temperature. It is noticeable that in both cases the 
two order parameters (77 and m) exhibit an anomaly at 
a given temperature, (a) T* ~ 0.70 and (b) T* ~ 0.97. 
These temperatures correspond to the bicritical and the 
triple point discussed in Figure I When crossing the 
triple point line (case (b)), the anomaly is accompanied 
by a discontinuity in the order parameters. 



IV. MONTE CARLO SIMULATION 



Monte Carlo simulations of model (10) have been per- 
formed in order to study the role of fluctuations. Start- 
ing from an initial (arbitrary) configuration, the subse- 
quent microscopic configurations are generated by using 
the standard Metropolis algorithm. We have focussed on 
two cases. First, on the stoichiometric alloy (c = 0.5) for 
different values of the parameters K* and J% ■ Secondly, 
we have fixed K* = 0.4 and J| = 0.6 and have studied 
the phase diagram as a function of c and T*. 



A. Simulation details 

The main results were obtained on a simple cubic lat- 
tice of size L=16(JV = LxLxL = 4096). More- 
over, a certain number of simulations with L — 24 and 
L = 32 were also carried out in order to study finite-size 
effects and to obtain illustrative real space snapshots of 
the system. Energy and order parameter fluctuations are 
measured according to the following definitions: 
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*2 



((H* 2 ) - (H*) 2 ) 



(19) 



(20) 



(21) 



The brackets stand for Monte Carlo (MC) averages, per- 
formed over a large number of uncorrelated configura- 
tions after the equilibration of the system. In order to 
find the phase diagram, the transition lines were located 
from the positions of the peaks of the above quantities. 
In many cases equilibration was checked by testing the 
fluctuation-dissipation theorem, i.e.: 



C = 



1 d{H*) 
N dT* ' 



(22) 



5 



Two kinds of numerical simulation experiments have 
been performed: 



the structure of the domains in the coexistence re- 
gion. 



Grand Canonical simulations. The simulations in 
the Grand Canonical Ensemble have the advantage 
of allowing faster equilibration. The alloy concen- 
tration is not fixed and an additional term taking 
into account the effect of the chemical potential dif- 
ference between both species is needed in this case. 
Formally, this is done by a Legendre transformation 
of the Hamiltonian (10). This yields: 



(23) 



A' 



Starting from an (arbitrary) initial configuration, 
the system at constant T* and fi* , evolves towards 
equilibrium by means of Glauber excitations pro- 
posed in both variables, tr^ and Si independently. 
The unit of time MCS (a Monte Carlo step) is 
defined as N independent proposals of each kind 
of flip on a randomly selected lattice site. Typi- 
cally the averages are performed over 1500 config- 
urations, taken every 20 MCS and discarding the 
initial 5000 MCS for equilibration. The regions of 
phase separation correspond to unreachable regions 
in the c — T* phase diagram. 

2. Canonical simulations. In these simulations the 
Glauber excitations are proposed in the magnetic 
variable Si only. In order to preserve the alloy com- 
position c, the variables (7$ evolve according to the 
Kawasaki exchange dynamics. The equilibration 
process is much slower in this case and the system 
may get trapped in metastable configurations. To 
get rid of such configurations, it is convenient to 
allow a certain fraction (q) of exchanges between 
n.n.n. atoms. Then, a MCS is in this case defined 
as N proposals of Si flips, N(l — q) proposals of n.n. 
exchanges and Nq proposals of n.n.n. exchanges. 
We have studied the effect of different values of q 
and found that q ~ 0.2 is enough to reach equi- 
librium in a reasonable time. Typically averages 
are performed over 3500 configurations, taken ev- 
ery 50 MCS, after discarding the first 25000 MCS 
for equilibration. In the region of phase separation 
the simulated system evolves to an inhomogeneous 
"slab" configuration with a flat interface. Because 
of finite-size effects, the energy of such configura- 
tions is very much dominated by the interfacial en- 
ergy and should be carefully analyzed. In spite of 
the long times needed to get reliable results, the 
simulations in the Canonical ensemble are very use- 
ful here since they provide information concerning 



B. Monte Carlo Results 

We start by presenting the transition temperatures as 
a function of the model parameters for the case of the 
stoichiometric alloy c = 0.5. This is shown in the lower 
part of Figure ^ (b) . A comparative look of both mean- 
field (a) and MC results (b) reveals that both solutions 
render the same qualitative behaviour. The fluctuations 
(taken into account in the Monte Carlo solution) have the 
effect of increasing the stability of the disordered, param- 
agnetic phases so that the overall transition temperatures 
are lower than in the mean-field solution. 

Figure [?] shows the /i* — T* section of the phase di- 
agram, drawn from the Grand Canonical simulations, 
with L = 16, J 2 * = 0.6 and K* = 0.4. We notice 
that the model parameters are those of Fig. |(d). It 
follows that both numerical simulations and mean-field 
techniques render the same qualitative phase diagram. 
The only remark that comes out is the smearing out of 
the re-entrant (OP) phase in the MC solution, due to the 
fluctuations. The available MC data does not allow for 
a conclusive determination of the nature (first or second- 
order) of the transitions. 

In order to compare data with experiment, the study 
of the c — T* section of the phase diagram is essential. 
It turns out to be a tough task because of the finite-size 
effects. In particular, to definitively resolve the coexis- 
tence region, one needs to use very large linear system 
sizes. 

Figure || shows the c — T* phase diagram correspond- 
ing to J| = 0.6 and K* — 0.4. In Fig. ||a we simulta- 
neously show the mean-field and the MC solutions. One 
observes that the main trends of both phase diagrams are 
the same. For practical reasons, we show the MC solution 
in more detail in Fig. ||b . The phase transition lines and 
the limits of the coexistence region, have been located 
from the peaks observed in the specific heat C. This cri- 
terion has been followed in both the Grand Canonical 
(open diamonds) and Canonical (black diamonds) sim- 
ulations . In the Grand Canonical simulations the co- 
existence region is revealed by unreachable zones in the 
(c) — T* diagram accompanied by flat steps in the curves 
of constant //* (three examples are depicted by small dots 
joined by a thin line). 

When comparing the results corresponding to the same 
system size obtained from both the Canonical and the 
Grand Canonical simulations we see that in the former 
the coexistence line occurs at lower temperatures. This 
is due to finite-size effects that have a strong influence 
on the stabilization of the mixed phase configurations. 
In this sense we have checked that when the system size 
is increased this effect is corrected and the phase sepa- 
ration occurs at higher temperatures. To illustrate this, 
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we have plotted in Figure g (with thick dashed lines) the 
upper part of the coexistence line obtained from Canoni- 
cal simulations, for two different values of the system size 
[L= 16 and L=24), as indicated. 

The same effect appears when studying the specific 
heat. In Figure |we show the temperature behaviour of 
the specific heat C (a) together with the order parame- 
ters to and r] (b) for J* = 0.6, K* = 0.4 and c = 0.25 
as obtained from the canonical MC simulations. Data 
shown correspond to L = 16 and L = 24. Note that 
the peak corresponding to the phase separation shows a 
much larger dependence on L than the peak correspond- 
ing to the order-disorder transition. The inset (c) shows 
the specific heat computed from the energy fluctuations 
(equation |l9|) and from the derivative of the average en- 
ergy (equation p2f ) . The agreement ensures that the equi- 
libration times considered are long enough. 

In spite of the difficulties described above, which cer- 
tainly hinder the location of the boundaries, the phase 
diagram presented in Figure H is essentially similar to 
that obtained experimentally (see Figure |l|) at least at 
moderate and low temperatures. The lack of resolution 
in the results makes it impossible to conclude whether or 
not the MC results render a line of triple points as oc- 
curs in the mean-field solution (lower part of [|). Unfor- 
tunately, the existing experimental data do not provide 
new information on this point. We suggest that more 
experiments are needed. Provided that the experimen- 
tal phase diagram is sufficiently well resolved, fine tuning 
of the parameters J| and K* (even Jj") would allow the 
matching of more details. 

Besides the determination of the phase diagram and 
the fluctuations, the MC simulations have the possibility 
of providing real space snapshots of the system configu- 
ration. Figure [l(] shows a two dimensional section of the 
simulated system with L — 24 for different homogeneous 
equilibrium phases corresponding to the phase diagram 
in Figure ||. Case (a) corresponds to the DP phase with 
c = 0.2 and T* = 1.2, (b) to the OP phase with c = 0.32 
and T* = 1.0 and (c) to the OF phase with c = 0.45 and 
T* = 0.79. The assignment of the different colours has 
been done by measuring the short-range order parame- 
ters in a cell of size 5x5x5 centered at each lattice 
site of a certain two-dimensional horizontal cut of the 
original system. When the values of the local magneti- 
zation m and/or local order parameter r\ are above 0.2 
the corresponding lattice site is considered to belong to 
a ferromagnetic and/or to an atomically ordered phase 
respectively. White, light gray and dark gray indicate 
DP, OP and OF regions. Black corresponds to local dis- 
ordered ferromagnetic regions which do not correspond 
to any stable phase. These appear because the fluctu- 
ations become both more probable and important with 
temperature in the homogeneous phases. Actually, the 
three snapshots correspond to a time evolution of 2 • 10 5 
MCS, when the average values of the long-range order 
parameters are perfectly equilibrated. Thus, the curved 
interfaces reveal that the fluctuations evolve with time 



and appear and disappear very quickly. 

In Figure [ll] we show snapshots of the system config- 
uration inside the coexistence region. The four pictures 
correspond to T* = 0.5 and to different values of the 
composition; (a) c = 0.1, (b) c = 0.2, (c) c = 0.3 and 
(d) c = 0.45. Note that for low concentration of the 
magnetic component (a), the OF phase consists of ferro- 
magnetic bubbles inside the DP matrix as expected. For 
larger values of c the ferromagnetic bubbles transform 
into rods or slabs (b). This is an artifact of finite-size 
effects that makes the system decrease the interfacial en- 
ergy by taking advantage of the periodic boundary con- 
ditions. Cases (c) and (d) are symmetric to (b) and (a) 
respectively. Given the large value of c, the matrix is 
ferromagnetic and the domains paramagnetic. 

It is known that the shape and the size of the mag- 
netic bubbles embedded into the non-magnetic matrix is 
crucial for the occurrence of magnetoresistance. In the 
light of the present results, we believe the present model 
is suitable for determining the optimum characteristics of 
such domains. Along these lines, the study of the kinetics 
of the domain growth after quenches from high temper- 
ature should supply useful information. This will be the 
subject of future work. 

V. CONCLUSION 

By using a simple lattice model we have shown that 
the magnetism of an ordered alloy may give rise to a 
low temperature phase separation between a ferromag- 
netic phase and a paramagnetic phase. The existence of 
this mixed phase is relevant in relation to the occurrence 
of phenomena such as paramagnetism and magnetoresis- 
tance. 

This study has been motivated by the behaviour ob- 
served in Cu3_ x AlMn x . Nevertheless, the strategy fol- 
lowed in the construction of the model should apply to 
other alloys. In particular, to those in which the fer- 
romagnetism is induced by the configurational ordering 
of the magnetic atoms, as occurs in Cua-^AlMn^. Our 
main conclusion is that this interplay between both kind 
of orderings is enough to produce the magnetic phase 
separation. We should mention that other effects such as 
elasticity due to the different atomic size of the elements 
may affect the final phase diagram. In spite of this and 
in view of the present results it is clear that the model 
captures the essential ingredients and makes it an appro- 
priate starting point for future dynamical studies of the 
kinetics of formation of the mixed phase after a suitable 
thermal quench. 
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TABLE I. Occupation probabilities for the stoichiometric 
DOz structure of CU3AI. 
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TABLE II. Occupation probabilities for the stoichiometric 
L2i structure of Cu 2 AlMn . 
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TABLE III. Guessed occupation probabilities for the 
non-stoichiometric DO3 structure with x > composition. 
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TABLE IV. Simplest occupation probabilities for the 
non-stoichiometric L2i structure with x < 1 composition. 
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FIG. 1. Approximate experimental phase diagram of 
Cua-xAIMn* from Ref. 24, 25 and 26 
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TABLE V. Guessed occupation probabilities for the 
non-stoichiometric C\iz- X AlMn^ alloys that include the DO3 
structure ( A = ) and the L2i structure (A 7^ 0) 
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FIG. 2. (a) Structure of the L2i phase of Cu 2 AlMn indi- 
cating the a, p, 7 and 5 sublattices. (b) Cell used for the 
present model. 
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FIG. 3. Mean-field fi*-T* diagrams for different values 
of model parameters JJ and K* , indicating the disor- 
dered-paramagnetic phase (DP), the ordered-paramagnetic 
phase (OP) and the ordered- ferromagnetic phase (OF). 
Dashed lines indicate first-order phase transitions while con- 
tinuous lines stand for continuous phase transitions. 
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FIG. 4. Mean-Field c-T* diagrams for two different sets of 
values of model parameters J| and K* , indicating the dis- 
ordered-paramagnetic phase (DP) , the ordered-paramagnetic 
phase (OP), the ordered-ferromagnetic phase (OF) and the 
coexistence regions. 
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FIG. 5. Behaviour of the order parameters r) and m with 
temperature for different values of c corresponding to the 
same two cases as in Fig. 4 
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FIG. 6. Dependence of the transition temperatures T* (□, 
■) and (°,») with J| for the stoichiometric compound 
(c = 0.5), from Mean-Field calculations (a) and Monte Carlo 
simulations (b). Open symbols correspond to K* = 0.4 and 
filled symbols to K* = 0.0. 
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J 2 =0.6, K=0.4, c=0.25 



FIG. 7. Phase diagram fi*-T* for J 2 * = 0.60 and K* = 0.40 
obtained from Monte Carlo simulations in the Grand Canon- 
ical ensemble. Lines are guides to the eye. 
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FIG. 9. Specific heat C and order parameters m and r\ as 
a function of temperature T* for = 0.6, K* = 0.4 and 
c = 0.25 obtained from canonical MC simulations. Data for 
L — 16 and L = 24 are shown with thin and thick lines respec- 
tively. The inset (c) shows the specific heat computed from 
fluctuations (continuous line) as well as from the derivative 
of the average energy (dashed line). 



FIG. 8. Phase diagram c-T* for J 2 * = 0.60 and K* = 0.40 
obtained from mean-field calculations (a) and Monte Carlo 
simulations (b). Thick lines indicate the phase boundaries 
between the homogeneous phases and the coexistence regions. 
The thin lines with dots correspond to Gran Canonical MC 
runs at constant chemical potential, diamonds indicate the 
positions of the specific heat peaks from Gran Canonical 
(open diamonds) and canonical (black diamonds) simulations. 
Thick dashed lines in (b) are estimations of the coexistence 
region boundary from Canonical simulations with L = 16 and 
L = 24. The continuous thin lines in (a) indicate the Monte 
Carlo phase boundaries for comparison. 
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FIG. 10. Snapshots of the system configuration cor- 
responding to the three different phases: (a) Disor- 
dered-Paramagnetic (DP) phase at T* = 1.2 and c = 0.2, (b) 
Ordered-Paramagnetic (OP) phase at T* = 1.0 and c = 0.32 
and (c) Ordered-Ferromagnetic (OF) phase at T* = 0.79 and 
c = 0.45. The shading identify locally the different phases 
according to the short-range order parameters as explained 
in the text. 
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FIG. 11. Snapshots of the system configuration in the co- 
existence region for T* = 0.5 and c = 0.1 (a), c = 0.2 (b), 
c = 0.3 (c) and c = 0.45 (d). 
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